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Abstract 

We consider relativistic U(l) gauge theories in 2 + 1 dimensions, with iV& species of complex bosons 
and Nj species of Dirac fermions at finite temperature. The quantum phase transition between the Higgs 
and Coulomb phases is described by a conformal field theory (CFT). At large iVj and Nf, but for arbitrary 
values of the ratio N^/Nf, we present computations of various critical exponents and universal amplitudes 
for these CFTs. We make contact with the different spin-liquids, charge-liquids and deconfined critical 
points of quantum magnets that these field theories describe. We compute physical observables that may be 
measured in experiments or numerical simulations of insulating and doped quantum magnets. 
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I. INTRODUCTION 



A number of experimental observations in the last two decades have begged for an understand- 
ing of interacting quantum systems that goes beyond the simple weakly interacting quasi-particle 
paradigm of solid-state physics. The large number of such systems with anomalous properties 
include as examples, cuprate superconductors, quantum critical heavy-fermion metals yj], Mott- 
insulating organic materials [2] and insulating frustrated magnets yfl. A new paradigm that has 
been invoked to describe these systems is quantum number fractionalization, see for e.g. Ref. 
Re-writing the original constituents in terms of fractional particles interacting with a gauge field, 
unconventional phases (or transitions) are accessed by the deconfined phases of these gauge theo- 
ries. These deconfined phases cannot find a consistent description in terms of quasi-particles that 
are compositionally related to the particles of the microscopic description, and hence provide a 
new paradigm that dramatically departs from the conventional confines of solid-state physics. 

An important issue that arises immediately is the detection of such exotic quantum phases or 
phase transition in experiments or numerical simulations. In most cases the appearance of these ex- 
otic field theories is also signaled by dramatic qualitative physical effects, e.g. a direct continuous 



transition between quantum states that break distinct symmetries [5] (such a transition would not 
be permitted in the conventional Landau theory). While these novel physical effects themselves 
signal the appearance of the deconfined field theory, it is desirable to have a direct understanding 
and numerical estimate for universal quantities associated with the continuum field theories. An 
estimate for the universal numbers and scaling functions allows for a quantitative comparison be- 
tween experiment/simulations and theory, and could be helpful for instance to distinguish between 
a non-universal weak first-order transition and a universal continuous one, in finite-size numerical 



computations [13 



14, 



15] 



The starting point of our analysis will be a relativistic continuum field theory of Nf flavors of 
charged fermions and N b flavors of charged bosons mutually interacting through a gauge field in 
two spatial dimensions. This continuum description makes no reference to the variety of micro- 
scopic lattice models that may realize the quantum field theory. Examples of the derivations of 
relativistic field theoretic descriptions starting from lattice models of spins or electrons in different 
contexts may be found in Refs. |^, [kI 11- We postpone a discussion of the application of our 
results to these physical realizations to Section PVT. 

We shall be interested in the general case of interacting field theories that involve bosons (z a ), 
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Dirac fermions (ty a ) and gauge fields (A^). Working in imaginary time the partition function of 
these theories may be written as functional integrals over the fields, 



Z = J V^ a Vz a VA^ex P {-S bf ) , 



(1) 



where the action S b f = Sb + Sf is specified below as the sum of bosonic and fermionic contribu- 
tations. 

In the field theory the bosons are described by the complex fields, z a , where a takes one of N b 
values. The iV& fields couple through their mutual interaction with a 'photon' gauge field, A^, in 
the usual minimal coupling that preserves gauge invariance. 

S b = -^ J d 2 rdT\(d^iA^)z a \ 2 (2) 

The constraint J2 a \ z a\ 2 = 1 must be satisfied at all points in time and space. The action Sb 
by itself is the well-known CP Wb_1 model. We remind the reader that in this theory, two-point 
correlation functions of z a itself are 'gauge-dependent' quantities, but correlation functions of 
the composite operator 4> a = z* a T^pZp are not. Indeed, the 'deconfined' theory of the Neel-VBS 
transition is described by Sb at Nb = 2, and 4> a plays the role of the Neel order parameter where 
T£g are the generators of SU(2), a physical quantity that mustbe gauge independent. Some large- 



181. See Ref.E 



for a review. We 



N b computations on this model are available in the literature H17L 
will recover these results as special cases of the more general field theory of interest here. 

We now turn to the action that describes the Dirac fermions. The fermionic part is described by 
the usual Dirac fermion action, 

S f = J d 2 rdr^ a ^(d, - iAJ]* a (3) 

where 7 M satisfy the Clifford algebra in 2 + 1 dimensions, {7^, 7 1 '} = 2S^ U . We use the isotropic 
version of Sf, since anisotropics have been shown to be irrelevant at the critical fixed point studied 



here 



fl. 



We shall consider the action Sb and Sf separately and also S b f = Sb + Sf in the large- 
Nb, Nf limit. In general, a Maxwellian term for the A^. Sa = J d 2 rdT(e afll/ d IM A u ) 2 should also be 
included. We note however that this term will not play an important role in our considerations: As 
we will show below, the contribution to the A^ propagator from the interaction with z a and ty a at 
large- N, decays with a slow power at long distances and hence dominate the contribution from Sa 
at large length scales. As is usual in quantum field theory, we introduce a finite-T by restricting 
the imaginary time integrals in the actions above to extend from to 1/T. The z a and A^ must 
then satisfy periodic and the ^ a anti-periodic boundary conditions. 



We note here that the action Sf, when considered separately is always quantum critical, i.e. 
there are no relevant perturbations consistent with its symmetry. This is in contrast to Sb which 
goes critical only at a single value of the coupling g = g c , i.e. the boson mass is a relevant 
perturbation at the phase transition between the Higgs (g < g c ) and the deconfined Coulomb 
phase (g > g c ). The results in this paper include the results for Sf alone by simply setting N b = 0. 
Results for the CP^ 6-1 model at criticality are obtained by setting Nf = and the general result 
for finite Nb and Nf applies for the Higgs to Coulomb phase transition in the background of critical 
Dirac fermion excitations. We assume that is non-compact throughout this paper so that the 
2+1 -dimensional Coulomb phase is unequivocally stable (by definition) to monopole proliferation 
(see e.g. Ref. ll6|) . 

The paper is organized in the following way: In Sec. |ll]we first introduce the formal machinery 
of the large-iV limit of Sbf considered here (large- iV in this paper implies large- Nb, Nf but an 
arbitrary finite value of the ratio). In Sec. Hill we then turn to an evaluation of the critical exponents 
and scaling dimensions of a number of quantities of physical interest at the T = critical point that 
separates the Higgs from the deconfined phase in the background of massless Dirac fermions. At 
the critical point of Sbf, certain susceptibilities (corresponding to rotations in the boson or fermion 
flavor space) have simple linear or quadratic temperature dependencies with universal amplitudes 
because of conservation laws. In Sec. [TV] we evaluate the numerical value of these amplitudes in 
the large- N limit. Finally, in Sec. |V]we provide a general discussion of the physical situations in 
which the quantities computed in this publication may be measured in experiments or simulations 
of insulating and doped quantum anti-ferromagnets. 

II. EFFECTIVE ACTION AND PROPOGATORS AT LARGE R, N f 

In this section we provide the framework of the large- Nb, Nf method used here and compute 
the effective action at large- N. At Nb, Nf = oo, fluctuations of the gauge field are completely 
suppressed and we have a theory of free bosons and fermions. Gauge fluctuation at next-to-leading 
order are conveniently computed by perturbation theory, since each photon line brings a factor of 
1/N. We hence require a computation of the photon propogator at order 1/N, in the theory Sbf 
introduced above at finite-T. In this paper we use two different gauges, depending on whether we 
are computing at T = or T ^ 0. At T = 0, when space and time are completely interchangable, 
it is convenient to use a gauge that reflects this symmetry, we use k^A^ fixed to a constant. At 
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(a) £V (b) D\ 



(c) Gq, (d) G z 




FIG. 1: Definition of diagrammatic symbols that arise from the action + Sf. (a), (b), (c) and (d), 
illustrate the four propogators that are used to construct the various graphs used here. See Eq. ([8]) for the 
T = values and Eqs. (1111121) for the T ^ 0. (e), (f), (g) and (h) show the four vertices allowed by our 
theory and their corresponding amplitudes. 

T ^ 0, for practical computation reasons, we prefer to use the gauge kiAi = (with % restricted 
to the two space indices). Throughout this paper, we will use greek symbols a, (3, fi, v etc. to take 
values x, y and r (2 + 1 dimensional space-time), and lower-case letters i,j, k to take values x, y 
(2-dimensional space). Summation is implied on repeated indices. 



A. Zero Temperature 

We begin with the resolution of the constraint on z a , by introducing a real field A, which acts 
as a Lagrange multiplier at each point of space and time. Including this into the bosonic part of 
the action, we can re- write it as, 

S b = - f d 2 rdr [Kfy - iAjzaf - *A(|z Q | 2 - 1)] (4) 

Now after a standard sequence of steps (see e.g. [16]) in the limit N b = oo the gauge field drops 
out and A takes on a uniform saddle point value that extremizes the action Ao = ir, 

d 3 p 1 1 



]% 3 p 2 + r g 
5 



(5) 



The propogator for the z particles is thus G z 1 (k) = k 2 + r. At ]V& = oo, it becomes critical when 
r = 0, i.e. at g = g c where j- = J 

The effective action is obtained to leading order in iV b , Nf by integrating out the bosons and 
fermions perturbatively, i.e. by evaluating diagrams shown in Fig. [2] At T = 0, it is possible to 
obtain simple closed-form expressions for the effective action, 

Sa-a = / ^A 2 + ^(V-^)^ (6) 



where 



U A (p,r) = N f ^ (7) 



+ N, 



p 2 + 4r / p \ y^r 

tan — = 

87rp \2y/r J An 



Details of this computation are provided in Appendix lAl We note that iT\(0, r) = {%ix^/r)~ l and 
11^(0, 0) = / -7, which is formally IR divergent in d = 3, but this will cancel out of all physical 
observables. 

Using the gauge k^A^ = 1 — ( for this T = calculation, we find the following form for the 
propogators, 



IV * q 2 

Dx = (XX) = 4" 

G 2 = (z*z) 



k 2 + r 

I 
k 2 



- — (8) 



Where we have used the usual notation that = fc M 7 M . The parameter ( is an arbitrary gauge- 
fixing parameter and it should drop out of any expression for a physical observable. 



B. Finite Temperatures 

At T 7^ 0, there is no particular advantage in using the conventional relativistically invariant 
notation of Eq. ©. So we use the following equivalent form of the fermion action, designed for 
transparent appearance of frequency sums 



S f = J d 2 xdr[¥ a (d T - iA T - ia x {d x - %A X ) - ia y (d y - iA y ))^ t 



(9) 




FIG. 2: One-loop diagrams that determine the effective action at large-iV. (a), (b) and (c) contribute to the 
effective action for the and (d) contributes to the action of the A field. 

Again, the index a — 1 ... Nf labels the fermion flavors, \I/ Q is a two-component Dirac spinor (for 
each a), a x ' y are Pauli matrices acting on the Dirac space. 

The general form of the effective action for the photon at large- A^, Nf can be simply written 
down by invoking a Ward identity, 



T 



d 2 k 



9 / 4 f 2 



(kiA 7 



tnAi) — h Aj Aj ( t\, , 



k 2 



D 2 (k,e r 



(10) 



where Di and D 2 are functions that can be evaluated at large- N by perturbatively integrating out 
both the fermions and the bosons. This process is illustrated in the Feynman-graphs in Fig. [2l A 
complete analytic evaluation at finite-T of Di 2 is not possible. However after some analytic ma- 
nipulations detailed in Appendix[A] we can bring the expressions into forms which allow efficient 
numerical evaluation of these functions, these evaluations will play an important role later in our 
T 7^ computations. In terms of Di 2 , we can immediately write down the photon propogator in 
the Coulomb gauge kiAi = 0. After imposing the gauge condition, the non-zero elements of the 
propagator are, 

1 



£>oo (q, e n ) 
Aj(q, e n ) 



Di(q, e n y 

m 3 



(ii) 



■i.i 



l 



q 2 ; D 2 ( q ,e n ) + (e 2 /q 2 )D 1 (q,e r . 
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The other propogators at finite-T are, 



^ z + k 2 + r 

G>$ = p (12) 

— %u n + a ■ k 

where k is the the vector (k x , k y ) and v m = l-nTm and tu n = 2nT(n + |) are the usual bosonic 
and fermionic Matsubara frequencies. 



III. CRITICAL EXPONENTS AND SCALING DIMENSIONS AT THE QCP OF S b / 

In this section we will work at T = and hence use the propogators derived in Sec. Ill Al 
and work in the gauge introduced there. As the coupling g is tuned, a critical point at which the 
condensation of z a takes place is crossed, for small g the system enters a phase with condensed 
z a and for large g the system is described by a 'liquid-like' state that breaks no symmetries (here 
the bosons z a acquire a mass and the low-energy theory is a strongly coupled 'algebraic' state 
of Dirac fermions coupled to a gauge field). In this section we first compute two independent 
critical exponents associated with the singularities at this transition (u N and r/ N ), including one- 
loop corrections about the Nf, N b = oo saddle point. We then turn to a calculation of the scaling 
dimension of the physical electron operator at the critical point. 

A. Computation of in S b f 

We begin with the correlation length exponent, u^, 

£ N oc (g - g c y UN , (defines v N ) (13) 

where £ N is the correlation length associated with the Neel order parameter. We note that in a 
given gauge one may calculate well-defined critical exponents using z a as the order parameter. 
The critical indices so calculated are in general gauge-dependent. We note however that since the 
correlation length is gauge invariant so is the exponent v z . Hence, un may be calculated without 
worrying about composite operators, i.e. = v z . Other exponents must however be calculated 
using the composite Neel field, we will present calculations for t]n in the following sub-section. 

We calculate the exponent of interest, u z , by calculating two gauge-dependent indices r\ z and 
7 Z and then using the scaling relation ^ z = v z {2 — r] z ). Since v z should be gauge-independant, any 
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(a) (b) 

. f7\ , . 



(c) (d) 




(e) - ® 



© €) 



FIG. 3: Self-energy diagrams that enter the two-point z a -boson correlation function at order 1/N. These 
diagrams are evaluated in a specific gauge to extract the gauge-dependent indices rj z and j z . A gauge 
independent index v can then be extracted from these quantities by using the exponent relation, j z = 

f«(2 - T) z ). 

dependence on the arbitrary gauge parameter, ( should drop out of the result. This is an important 
check on our results. We begin by defining the anomalous dimension of the z a field: 

dim[;2 Q ] = 2~~~~^~' ( m a gi ven g au g e ) defines r] z ) (14) 

We will make the gauge-dependence of r/ z explicit by calculating its value for arbitrary (. The 
calculation of r] z is straightforward and is obtained by picking up the co-efficient of the p 2 logp 
pieces in a perturbative expansion of G z , and then re-exponentiating. Diagrams in Fig. (3^a,c) 
contribute giving, 

Th ^. 7Tr A T7 ^\^ +2( \ (15) 



37T 2 N b (N f + N b )7T 



17. 



3 S 



For ( = 1 and Nf = 0, this reproduces the result of Ref. 

The critical index ^ z determines how the mass of z a (in a given gauge) goes to zero at the 
critical point, 

G-\k = 0) = {g- gc) y % (defines lz ) (16) 
To leading order we may just use r c = E(0, 0) and ^ = J we find: 



Introduce a new parameter r g , 



1 1 



9c 9 J P P 2 P 2 + r g 4vr 
which can be related to r by, 



(18) 



4l £"/>>.<» C9) 

Beyond the ]V& = oo approximation the critical point is no more at r c = 0. In general, 

G-\p) = (p 2 + r) -E(p,r) (21) 

and so, r c = E(0, r c ). 

Now -writing the boson mass in terms of the new parameter r g , 

G-\0) = r g - (s(0,r g ) - ^Me(0,0)) , (22) 

we have eliminated the quantity r completely. The goal now is to compute the term in brackets by 
evaluating the self-energy diagrams and pick up the terms with r g log(r ff ) divergences. Then these 
logs must be re-exponentiated yielding the exponent j z : 



i-ii 



^1 + a log 


"A 2 " 






^ ~ \g-9c\ 




J 9. 



|2(l-o) 



(23) 



where in the last equality we have used the fact that r g oc \g — g c \ 2 , as derived above and that a 
vanishes in the iV — > oo limit. 

In order to demonstrate the usefulness of these formal manipulations, we calculate the index 7 A 
for the 0(N = 2N b ) cr-model. This index is of course completely gauge invariant and follows by 
including only the A contributions to the self-energy, i.e. the diagrams in Figs.[3£c,e). 

E ( C ) : 2* 2 f d 3 q 1 1 



2! J 87r 3 n A (g,r) (p + q) 2 + r 
M _ 41i< 1 fd\ 1 

4! n A (0,r) 7 87r 3 (g 2 + r) 2 1 ; 

d3p 1 1 (25) 



87r 3 I1a(p, r) (p + q) 2 + r 

Now putting these expression back into Eq. d22l) . we find that ^|°'°| s(0, 0) = 0. Note however 
that this term cannot be ignored in general. Indeed, when we include gauge fluctuations below, we 
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find that it does not evaluate to zero. In the current case, because of the vanishing of this term, we 
find that the quantity in brackets in Eq. (1221) is simply £( c )+( e ) (0, r), 



d 3 p 



8vr3 n A (p,r) 



1 + n' A (p,r) 



j? + r 2n A (0,r)_ 



(26) 



where n' A = 5^-. Now expanding the integrand for large p we can extract the co-effecient a A 



■^Tf^- Thus 7 = 2 — -^f, reproducing the standard result for the 0(N) model 112 IN . 

We now turn to the inclusion of the gauge field: There are four diagrams [Fig. [3ta,b,d,f)] that 
contribute to = E^) + E^ + E^ + E^ . Including all factors of i 2 , contractions and factorials, 
they are: 



E («0 



d 3 q 



1 



iV(p, g) 



2! J 87r 3 n A (g,r) (p + g) 2 + r 

( " 3 + C) /0n^7) 

4i 2 1 f d 3 q 1 
2!2!n A (0,r)/ 8^n A (g,r) 



x 



E (/) 



^7r 3 (p 2 + r) 2 ((p + g) 2 + r) 
(-3 + C)2i 2 1 f d 3 p I 



x 



2! 1! n A (0,r) J 8vr 3 (p 2 + r) 2 
d 3 q 1 



(27) 



8tt3 n A (g,r) 

where iV(p, q) = 4p 2 — ^(p • g) 2 + (4p • g + g 2 )(l — £), f° r a general gauge. Note that generally, 
£(&) + £(/) = o, so they never make an appearance. 

Our task is now to compute Eq. ((22)) with E = E^ + E^ . This expression seems to be plagued 
with both IR and UV divergences, we will see below that both divergences exactly cancel as they 
must. 



E(0,r) 



q Ii A iq, 0)q 2 + r 
1 



;i-o 

N(p, q) 



1 



E(0,0) 



Il A (0,r) J pq (p 2 + r) 2 [(p + q) 2 + r] U A (q,r) 

' (i-C) 



ru(g, o) 



l 



N{p, q) 



1 



n A (0,0) L qP 4 ip + q) 2 Tl A iq,0) 



(28) 
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Using these in Eq. (O, we find for the term in brackets, T ( ... } = S(0, r g ) - ^^S(0, 0) 



(...) 



+ 



(i-C) 
1 

n A (o,r,) 
1 

n A (0,r 9 ) 



U A (q,r g ) q 2 + r g 

N(p, q) 



(29) 



p,q 



(p 2 + r g ) 2 [(p + q) 2 + r g ] Tl A (q, r g 



N{p, q) 



;i-o 



P A n A ( g , 0) 



Xp + q) 2 

Note that the last term was formed by combining two terms, both of which were IR divergent, this 
term however has no IR divergence. No other term has a IR divergence. We now proceed to show 
that all the UV divergences also exactly cancel. In order to do so, let us re- write, 



(...) 



;i-o 
i 



^A(q,r g ] 

1 



U x (0,r g ) J q U A (q,r 



+ (i-C) 
— J A(g,r 9 ; 



H-A(q,r g ) q 2 + 



+ 



1 



n A (0,r,)7 4qU A (q,0) 



(30) 



The following integral was used on the last line: J p 
Now, we can rewrite this integral as, 



N(p,g) 

(p+q) 2 



:i-o 



T<... 



(...) 



q 2 dq 
1 



C 



H-A(q,r g )q 2 + 



(3D 



lA(q,r g 



+ 



1 



1 



II A (0, r g ) Il A (q, r g ) n A (0, r g ) Aq U A (q, 0)J 

All we have to do is to expand the term in [...] for large q to extract the log divergence (Evaluation 
of I A is detailed in the AppendixEl There should be a 0(l/g 3 ) power but all lower powers must 
cancel. This is indeed found to be the case: 



T, 



( ...) 



q 2 dq 



2tt 2 
a A r g log 



7N f -W b 



c 



g (N b + N f ) 2 ' N f + N b 



1 



A 2 



(32) 



where a A = -^( ^^ + 

Using a = a\ + a A , and 7 = 2(1 — a), we have the final expression for 7: 



7 



12 8 f7N f -9N h 
+ 1 



7T^ 



+ 



c 



N b 7i 2 \(N b + N f ) 2 N f + N b 



(33) 
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We note that when £ = 1, Nf = 0, this reduces to the result of Ref. 17 (they only present results 
in the Landau gauge, £ = 1). Using their quoted values of 77 and v and 7 = u(2 — i]), we find that 
their result gives (one has to also adjust N by a factor of 2, because they have a theory with N/2 
complex fields): 7 = 2 — 76/(ir 2 N b ), in agreement with us. 

We can now calculate the co-efficient v by using a scaling relation: 



Iz 



2 K 2' 



(34) 



1 - 



TT 2 N 



+ 



4 f7N f -9N, 



7T- 2 



+ 



c 



1 + 



.10 



37r 2 N b (N f + N b )ir 2 1 3 



1 - 



16 



+ 



4 7N f -9N b 



+ 2C] 



20 



3n 2 N b tt 2 (N b + N f ) 2 (N f + N b )3ir 2 
Reassuringly, the gauge parameter £ drops out as expected. Also for Nf 



consistent with Ref. 



17 



0, v 



48 



B. Evaluation of i]n in S b f 

We now turn to an evaluation of the scaling dimension of the Neel field. We first define the 
anomalous dimension r/ N through the equation: 

a- u«i D-2 + r] N 

dim^] = (35) 

where (p a N = z* a T^z p is the Neel field and T a are the generators of SU(iV). Note that at N f ,N b = 
00, the index r]N = 1, i.e. it is already non-zero. This is unlike, e.g. the O(N) model, which at 
N = 00 has no anomalous dimension. Since the Neel field is a composite of two z a fields, its 
scaling dimension may be written as, 

dim[0^] = 2dim[z a ] + 77 vrtx (36) 

where i] vvtx gets contributions exclusively from vertex contributions illustrated in Fig.|4|a,b). This 
contribution is conveniently calculated by including a source h a cj) a N in the action and studying its 
renormalization, by momentum shell RG. There are two contributions coming from the photon 
and A propogators (integrals are from A/s to A). They give: 



2 J (2vr) 3 n A g 4 2 J (2vr) 3 1I4 g 4 



7T 2 N b &W 7T 2 (Nf + N b ) 
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FIG. 4: Vertex corrections to the scaling dimensions for the composite operators: (a,b) Neel vector, <p a N 
(c) electron operator, c = z*^ and (d) fermion bilinear 



This gives ?7 vrtx = - ^(^+^ (1 - 0- We thus have the result, 



r] N = D - 2 + 2r] z + 2r] VItx 
32 128 

= 1 + 



3n 2 N b 3n 2 (N f + N b ) 
Again ( drops out as expected, since rj N should be gauge invariant. 



(38) 
(39) 



C. Scaling Dimension of the Electron Operator at Criticality 

The scaling dimension of the physical electron operator is calculated at large- N by identifying 
its scaling dimension with the bilinear c = ^!z. The scaling dimension of c must clearly by 
gauge-invariant. 

The scaling dimension of * is defined by: 

dim[»] = B^JhL (40) 

We can evaluate 77* by evaluating the lowest self-energy graph S(fc) due to the gauge field and 
then picking up the )k\og(k) part: 



2! J (2vr) 3 n A '^(k + q) 2 \ ^ * q 2 

1 ^log(fc) 



3n 2 (N f + N b ) n 2 (N f + N b ) 
This gives the gauge-dependent result, 77* = 3n 2 (N 8 b+Nf) - n 2 { * +Nf y 
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Now we can get the gauge invariant quantity, dim[c], by including the vertex contribution in 
Fig. He). 

dim[c] = dimfz] + dim[*] + r? vrtx (41) 

(42) 

The parameter £ drops out as expected. Evaluating the vertex contribution by the momentum shell 
method, we find, 

l-C)log(s) (43) 



7T 2 (JV> + N b 
giving the final answer, 

dim[G] = 2 + 3^ " 3n*(N b + N f ) ^ 
Finally we note that it is possible to construct yet another gauge invariant bilinear, ty^ty. Its 
scaling dimension, dim[\[>t\l/], can again calculated for arbitrary £, we find 

64 

dim =2 (45) 

L J 3<K 2 (N b + N f ) 

This quantity too gets a vertex contribution, Fig. Htd), in addition to the 2dim[\P] part. Again the 



full result is reassuringly gauge-invariant and for N b = reproduces the results of Ref. 
calculated this number in the context of the 'algebraic spin liquid'. 
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IV. UNIVERSAL AMPLITUDES OF SUSCEPTIBILITIES AND SPECIFIC HEAT AT FINITE-T 
FROM FREE ENERGY 

In the previous section we studied the action S bf at T = 0, which corresponded to an infinite 
3-dimensional system. We now turn to a study of some properties at T ^ 0, which corresponds 
to studying the field theory in a slab geometry in which two dimensions are infinite, but the third 
direction has a finite extent of 1 /T. This naturally complicates the computations, since the rela- 
tivistic invariance of T = is destroyed and frequency integrals are replaced by discrete sums on 
Matsubara modes. In this section we evaluate the temperature dependence of the specific heat and 
certain susceptibilities that correspond to rotations among the fermionic and bosonic fields respec- 
tively. All quantities discussed in this section may be computed exactly from the low-temperature, 
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low-magnetic field (to be defined below) expansion of the free-energy. We compute such an ex- 
pansion in a the large N b and Nf limit, allowing for an arbitrary value of the ratio, Nb/Nf. Since 
the computations are done at finite-T we will work in the gauge introduced in Sec. Ill Bl and use 
the form of the propogators derived there. 

We shall consider two imaginary "magnetic fields", H b and Hj that are applied in the "di- 
rection" parallel to the SU(N) generator, diag(l 1 1 ... -1 -1 -1...). The chief advantage of us- 
ing this generator is that it modifies the action in a very simple way. For the bosons, it shifts 
uj n — > u n — 6H b /2, where 9 = 1 for the first half of the components of z a and 9 = —1 for the 
others. For the fermions it does the same, with uj n — > oj n — 9Hf. 

Once we have the free-energy as a function of H b , Hf and T, we can extract the specific heat 
and susceptibilities by differentiating, 

d 2 T 

Xb = M! (46) 

d 2 T 

Xf = jnn (47) 



OH 2 

C - 9 



dTT 



d 2 T 

■T m <«> 



n 

From the analysis of Ref. |22L it is know that Xb, Xf ^ T and Cy oc T 2 . These simple integer 
scaling powers are due to energy conservation and bosonic and fermioinic rotational symmetry. 
The proportionality constants of CV,x&, and Xf are universal amplitudes (^4.^)^X6 an( ^ ^x/)' 
divided by the square of a non-universal velocity, c. Since c is non-universal, it must be determined 
separately for each model. In the field theoretic analysis of this Section, we will set c = 1, though 
its presence must be kept in mind for comparisons to simulations/experiments. 

At large and Nf we can expand the free energy, T = —T In Z 

T = N b f b + N f ff + + f A (^j . (49) 

Note that in the N = oo limit, the gauge field disappears and the bosons and fermions de-couple 
so they contribute separately; in next to leading order this is not true. Hence, f ob and f°f are 
independant of the ratio N b /Nf, whereas f 1A is a function of the ratio (f ix can only depend on N b 
and hence is also independent of the ratio N b /Nf). In order to calculate the universal amplitudes 
associated with the specific heat and the susceptibilities we compute f ob , f°f, f ix and f 1A for 
arbitrary T,H b ,Hf and N b /Nf. We first outline the computation of f ob ,f°f and f 1 , then we 
present the numerical results for each of Xb, Xf and C v , as the ratio N b /Nf is varied. 
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When N b , Nf = oo, we only need to consider a simple gaussian theory of bosons and fermions 
and the computation of the free energy is straightforward. Including the shift in uo n as discussed 
above, we obtain for f ob and f°f, 



m 
9 



f ob = \ E / H e + ( w » + ^/ 2 ) 2 + ™ 2 ) + ln ^ 2 + ( w » - ^/ 2 ) 2 + m2 )] " 

/°' = -T y [ln(l + e -( fc+J ^)/ T ) + ln(l + e -^ H ^ T )] (50) 

where the mass, m of the bosons which at the critical coupling depends only on T and the applied 
field H b . It can be computed by minimizing the free energy at N b = oo. This leads to the following 
non-linear equation. 



L / 4^2 



1 



1 



2 k 2 + (oo n + 6H b /2f + m 2 



rf 3 p 1 

87T 3 p 2 



(51) 



where ^ should be summed over ±1. Applying the Poisson summation formula and evaluating the 
integrals and sums, we have the equation for m, 



d 2 k 



1 



1 



e (Vk 2 +m 2 +ieH b /2)/T _ I 



m 
An 



which has the solution 

m = T In 



4tt 2 2y/k 2 + m 2 

2 cos(H b /2T) + 1 + v/(2cos(iJ 6 /2T) + l) 2 - 4 



(52) 



(53) 



Computation of the various quantities of interest at N = oo follows simply by differentiating the 
free energy. 

Turning now to the 1/N correction to the free energy, fi receives two contributions, one from 
the Gaussian fluctuations of the A field and the other from the Gaussian fluctuations of A^. The 



result is, 



/ 



1A 



° 2 + t 2 ° l 



(54) 
(55) 



where n A , Di and D 2 axe evaluated at (k, e n ), and it must be kept in mind that they depend on 
H b: Hf through the shift of all internal boson and fermion frequencies (as detailed above) in the 
evaluation of the one loop graphs for the A and propogators. We note that the value of g c , gets 
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N b /N f 

FIG. 5: Gauge field fluctuation contribution to the specific heat amplitude, A}f plotted as a function 
Nb/Nf. The values at Nb/Nf = 0, oo are shown for reference as dashed lines. 

a correction from its N = oo value, j- = j (see Sec. Ill A)) , which can be evaluated from 
Eq. (|2TT) to be Sj- = N + Nb J ;p> tms shift which is cut-off dependent must be included in 
the free-energy to properly cancel all the UV diveregences at 0(1/ N) and produce the universal 
cutoff-independent results described below. 

We now turn to a presentation of our results for CV, Xt and Xf> obtained by numerically differ- 
entiating computations of the free energy formulae outlined above. 



A. Specific Heat 

The quantum critical specific heat at low-T is of the form C v = Ac v T 2 , where we can organize 
the large- N expansion of the universal amplitude, Ac v as, 

A Cv = N b A° c b v + N f A° c f v + A l c x v + A l c A v (N b /N f ) (56) 

Differentiating Eq. (gO]), we obtain the estimates, A°c = 6^ w 1.83661 and A°J = 6^) 



1.72182. Note that the proportionality of Aq to £(3) with a rational co-efficient is non-trivial [19]. 

Turning to the 1/N corrections, the computation of A}$ and A}^ are far more com- 
plicated and involve some tedious numerical calculations. Basically numerical computations 
of the finite temperature propogators (as detailed in Appendix |A]) have to be inserted into 
Eq. (l54l) which is then itself evaluated numerically. From our numerical analysis, we obtain, 
| J2 n f In (8^fc 2 + e2II(fc, e n )) = -0.03171T 3 , allowing us to estimate, Ag v = -0.38368. 
We can get A]^. from a similar analysis, we plot it as a function of N b /Nf in Fig. |5J 
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FIG. 6: Gauge field contribution to the amplitude of the susceptibility to applied H b as a function of Nj /N b . 



B. Susceptibility to applied H b 



The susceptibility to H b is obtained by differentiating the free energy twice with respect to 
Hb (keeping Hf = throughout). The large— N expansion of the amplitude A Xb defined from 
Xb = -A. Xb T can be organized as, 



A Xb = N b Al + A x x b + A x A b {N f /N { 



(57) 



We obtain the estimate A° x b b = In ( j w 0.0856271. From our numerical analysis described 
briefly above, we find f En / fl ln + e^fc, e n )) = -0.03171T 3 - 0.0132577 fe 2 T. This 

allows us to estimate w —0.02650. Finally from the numerical evaluation of the gauge field, 
we plot the function A l x A (N f /N b ) in Fig. |6] 



C. Susceptibility to applied Hf 

Finally, we turn to an evaluation of the susceptibility to Hf which also has a linear-T depen- 
dence. The universal amplitude appearing (just as in the case of Xb above) can be expanded as, 

A Xf = N f A% + A x A (N b /N f ) (58) 

In this case we obtain A° x f f = 2^ w 0.220636. A computation of the function A lA {N h /N f ) 
along the lines described above is plotted in Fig. [7] 
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FIG. 7: Gauge field contribution to the amplitude of the susceptibility to applied Hj as a function of 

N b /Nf. 

V. CONCLUSIONS 



All the computations carried out in this paper have been for the continuum field theory, and 
are hence of interest to any problem which realizes this field theoretic description. We now list 
three different physical problems of interest to which these results are applicable, and provide a 
brief summary of the connections between physical quantities and the universal numbers we have 
computed. 

Deconfined Neel-VBS transition: The deconfined Neel-VBS transition has been shown to be 
described by the field theory Sm with Nf = 0. The results of Sec. [Ill] apply directly to the 
critical exponents of this transition (which has Nf = 0, Nb = 2). In particular, Eqs. (|34l) and (1381) 
characterize the two-point Neel correlation function. At the Neel-VBS quantum critical point, 
the uniform magnetic susceptibility, Xu should depend linearly on temperature. This quantity is 
exactly the Xb, Eq. (1571) computed in Sec. [IV] with Nf = 0. The result for Cy, Eq. (1561 ) with 
Nf = also applies to the deconfined Neel-VBS transition. 



Algebraic Charge Liquids: In a recent publication ft 1 111 , a new class of states of matter, 



dubbed 'algebraic charge liquids' (ACL) were shown to arise naturally in doped quantum anti- 
ferromagnets, which are close to the deconfined Neel-VBS quantum critical point. In the language 
of this study, z a are bosons that carry spin and ^ are the superconducting quasi-particles. In this 
study the transition from a superconducting Neel state to the superconducting ACL state is de- 
scribed precisely by the field theory under study here with N b = 2 and Nf = 4, whereas the 
criticality of the ACL state is described simply by iV& = 0, Nf = 4. The critical exponents of 
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Sec. [in] apply directly to the transition from the superconducting Neel state to the s-ACL state. It 
has been shown [11], that the susceptibility Xf discussed in Sec. |IV]is related to the finite-T con- 
tribution to the superfluid density of the s-ACL state (the result for Xf at N b = was published in 
Ref. 11). Just like the Neel-VBS case the susceptibility Xb is the uniform magnetic susceptibility, 
and in Sec. [IV] we have the results for the universal amplitude both at the critical point and in the 
s-ACL phase. 

Algebraic Spin Liquids: Finally the result for Xf with Nh = applies to the fermionic spinon 
theories of algebraic spin liquids of square [23|, kagome and other lattices. In this case, Xf is 
related to the uniform magnetic spin susceptibility because the fermions carry spin. 
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APPENDIX A: EFFECTIVE ACTION AT LARGE-N: EVALUATION OF D x , D 2 , ILa AND n A 

In this section we provide the details on the evaluation of the functions D\ and D 2 that de- 
termine the effective action for the photons at large- N. We evaluate the general form of these 
functions at large N b , Nf and at finite-T, as defined in Eq. (flOl) . From these results we can get 11^ 
relevant to Sec. Ill Al by setting T = and completing all integral analytically. Finally we provide 
a computation of the A propogator. 

The evaluation of D\ and D 2 corresponds to evaluating the diagrams in Fig. |2]at finite-T. Each 
functions gets a contribution from the bosonic as well as fermionic loops. We evaluated these 
separately, writing D t = D lb + D\f and D 2 = D 2b + D 2 f. We first present an evaluation of the 
fermionic loops and then turn to a computation of the bosonic contributions. 
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1. D\f,D2f: Fermionic Loops 



The only diagram that needs to be evaluated for the fermions is the one appearing in Fig. [2tc). 
Starting with the function Dif(k, e n ): 

= -N f / Tr[G^(q,u; n )G^(k + q,^ n + e n )] 
uj n (uj n + e w ) - q ■ (q + k) 

> 2 + g 2 )(K + e n ) 2 + ( q + k) 2 ) K } 

-2 {2u n + e n f + k 2 

q 2 + + {ul + q 2 ){{u n + e n ) 2 + (q + k) 2 ). 

We will use the shorthand J = T J to signify a summation on the internal frequencies 
and momenta in this paper. Evaluating the first term by ^-function regularization, and performing 
the momentum integration in the second term using the usual Feynman trick, we find, 





£>i/(k, e n ) 



N f T In 2 : N f T 

7T ATT 



+ 



iV f T In 2 N f T 

— + -7— 

7r Arc 



E 

E 2 



dx- 



{2cu n + e n ) 2 + k 2 



1 - x)uj 2 n + x(u n + e n ) 2 + k 2 x(l - x) 

(2uj n + e n ) 2 + k 2 A + k 2 + uj 2 n + (w n + e r 
— m — 



.4 



2|cj„(u;„ + e r 



where 



A 2 = (k 2 + (|w n | - \u n + e n |) 2 )(k 2 + (\u n \ + |w n + e n |) 2 ). 



(A2) 



The frequency summation is formally divergent, but we can subtract the divergent piece by zeta- 
function regularization to yield 

N f T In 2 



D lf (\s.,e n ) 



+ 



7T 

Arc 



E 2 



(2w n + e n ) 2 + k 2 A + k 2 + u 2 n + (w n + e r 



.4 



In 



2|o; n (c t ; n + e n ) 



-4 



(A3) 



This frequency summation is now evaluated numerically. For this, it is useful to know the large cu n 
behavior of the term in the curly brackets. After symmetrizing the positive and negative frequen- 
cies we obtain: 



(k 2 + e 2 



3w 2 



+ 



94 - k 2 

30u;4 



50e 4 - 19e 2 k 2 + k 4 
210^1 



+ . . . 



(A4) 
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The asymptotic behavior is then summed by using the identities 

V 1 1 1 7 

^ (In - l) 2 ~ AN ~ A8N 3 + 960N 5 4 ' " ' 

n=N+l v ' 



5^ (2n - 1' 



1 

+ 



(2n - l) 4 A8N 3 96N 5 

n=N+l v 7 



oo 



(2n - l) 6 32(W 5 

n=N+l v y 

We now turn to a similar evaluation as above for D 2 f, 

,2 r Ai 



k 2 /" d Q 

D 2f {k,t n ) = --r-r NT z2 / ■^ Ti: [ 'xG^(c[,uj n )a y G^(k + q,u n + e n )] 
— 2NT^ / rf2g Zqxqy + q x k y + q y k x 

xk y u n * 

7T 7 ^ ((1 - x)u; 2 + x(u n + e n ) 2 + A; 2 x(l - x)) 



k x k v ^ J 4tt 2 {ul + q 2 )(K + e„) 2 + (q + k) 2 ) 



iVf T f , (cu ra + e n ) 2 - n> 2 (q; w + e n ) 2 

7T ^ I 2k 2 11 Ul 



n 



((o; n + e n ) 2 - ujjf + k 2 ((cu n + e w ) 2 + u; 2 ) 
Ak 2 

xln^ + ^ + ^ + ^ + ^U (A6) 

2|^„(^„ + e n )| J 

The large uj n behavior of the term in the curly brackets, after symmetrizing the positive and nega- 
tive frequencies, is: 



1 7e 2 - 2k 2 13e 4 - 34e 2 k 2 + 3k 4 

j iL j fl It j 



Qui 6(W 4 420^ 
At T = 0, the exact values of -Di,2/ can be computed analytically, 



(A7) 



D lf (k,u) = D 2f (k,u) = ^^==. (A8) 
2. Du, D 2b : Bosonic Loops 

We now present the evaluation of D lb and D 2b that arise from the diagrams in Fig.[2a,b). The 
function D lb is 
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D lb (k, e r 



N b Tj2j 



4tt 2 



q 2 + o> 2 + r 
(2u n + e n f 



(q 2 + o> 2 + r)((q + k) 2 + (w„ + e n ) 2 + r) 
The first term we evaluate by dimensional regularization 

d d q 2 



47T 2 q 2 + a; 2 + r 



(cu 2 +r) 1 - d / 2 

2-d 



T 



r(i-d/2) v 

( 4vr )d/2 2^ 

r(i-d/2) ^ 

( 47r )rf/2 

+ (d-2)X>(l + -L) 

n=l V n 7 



(27T) 



2-d 



C(2-d) 



1 

27 



K£) + |>( 1+ *) 



(A9) 



(A10) 



(All) 



(A12) 



where we have used the large-iV value for the mass parameter y/r = 2 m((\/5 + l)/2) in the last 



line (see Ch. 5 of Ref. 



20| for more details). 



D lb {k, e n ) 



47T 



E 



dx 



(2o> n + e r 



(1 - x)u)\ + x{uj n + e n f + r + k 2 x(l - x) 



(A13) 



The frequency summand has the following expansion at large o> n 

2k 2 + 12r - e 2 n -4k 4 + k 2 (-40r + I7e 2 n ) - 3(40r 2 - 50re;; + 3e 4 



+ 



3a; 2 30a; 4 
6k 6 + (84r - 73e^ ) k 4 + (420r 2 - 826e*r + 81e 4 ) k 2 + 10 (84r 3 - 315e 2 n r 2 + 147e 4 r - 5e£) 



210^ 



For Do?, we have 



^2fe(k, e n ) 



N b T 
An 



k 2 (2x- if 



o (1 - x)u 2 + x(u n + e n f + r + k 2 a;(l - x) 



(A14) 



and the following expansion at large cu n 



k 2 + lOr - lie* k 4 + 70r 2 - 266re;; + 86e 4 + 14k 2 r - 23k 2 el 



3a; 2 



30a; 4 



+ 



210o;6 



(A15) 
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3. U A and n A : Action at T = 



The above expressions for Di, D 2 at finite-T were brought into a form from which efficient 

numerical evaluation was possible. At T = 0, a full analytic evaluation of the integrals entering 

the effective action is possible. We will outline the steps of these calculations now. We begin with 

a computation of II A . Just as in the finite-T case, U A = Y\ b A + 11^ receives contribution from both 

bosonic and fermionic loops, which follow from the T — > limit of the above expressions for the 

finite-T action. Using the form in Eq. ©, we can deduce by looking for example at the q x q y term, 

_ u b<h% = _ N f d 3 q (2k + q) x (2k + q) y 

A q 2 b J 8tt 3 (k 2 + r)[(k + q) 2 + r] ^ 

Now the evaluation of this integral follows the standard steps (outlined for example in Appendix|B]) 
and we find, 



n b A (p,r) = N b 



r p 2 + 4r _-, / p \ Jr 
■ tan 



(A17) 



8-rcp \2Vr J 4 7r 

The fermionic contribution can be evaluated in a similar way or by noting the simple relationship, 

n A = E ^D 2 (p, co), giving from Eq.|A8]that, 



U a(p) = (A18) 

Putting these results together we get the expression for IT^ quoted in Eq. ©. It is curious to note 
that when r = the bosonic contribution becomes H b A (p, 0) = exactly the same as for the 
fermions. 

The evaluation of I1a follows the same steps as for except with different contributions at 
the vertices, Fig.[2td). 

/d 3 k 1 
a? (g + r) [(* + ,)» + ,•] (A19) 

which can be evaluated using the usual Feynman parameter and integration shift of Appendix [Bj 
resulting in the form quoted in Eq. ©. 



APPENDIX B: EVALUATION OF I A 

Here we explicitly show how to calculate I A that appears in computations of Subsection IIIIAI 
We note that we have used the standard sequence of steps illustrated in this example below many 
times in the course of the computations presented in this paper. 

d s p V - %P ■ qf + (4p • q + q 2 )(l - () 



57T 



3 



{p 2 + r g ) 2 [(p + q) 2 + r g 



(Bl) 
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FIG. 8: Diagrammatic contributions to the susceptibility, Xf arising (a) at N = oo, (b,c) at order 1/N. All 
the diagrams give a T-linear dependence as expected for the susceptibility of a conserved charge for z = 1 
criticality in d = 2. 

The integral is evaluated by executing the following steps in sequence: introduce a Feynman 
parameter x, shift to I = p + xq, integrate over / and then integrate over x. 



fc = jf*2<l-,)/g| 



d 3 p V - f(p • q? + (4p ■ q + q 2 )(l - () 



dx 2(1 — x) 



dx2(l - x) 



x)(p 2 + T g ) + x[(p + q) 2 + r g ]) 3 
dH 4/ 2 (l -§) + g 2 (l - C)(4s(z - 1) + 1) 

8^3 (/2 + r + q 2 x (l _ x ))3 

l-C/3 | [4x(x-l) + l](l-Qg 2 



^tt a/ ? 2 (1 — x)x + r 



32tt 



^A/g 2 (l — x)x + r 



1 



;i-C)n A (o,r 5 ) + - 



4g vrg^ 



+ 0( 3 ) 



9 



(B2) 



The neglected terms produce UV convergent diagrams when the integrals over q are evaluated and 
hence can be neglected. 



APPENDIX C: DIAGRAMATIC EVALUATION OF FERMION SUSCEPTIBILITY 

In this Appendix we calculate one of the universal amplitudes of Sec. [TV] A Xf , by an entirely 
different diagrammatic method. The perfect agreement between the two methods provides a non- 
trivial check on our rather involved computations. 
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Consider the susceptibility, x.f of the SU(iV/) charge, Q a = v]>t T^Jifp, where T a is a generator 
of the symmetry. Since the SU(iV/) charge is a conserved density it acquires no anomalous dimen- 
sion. Indeed a simple scaling analysis |22J] establishes that for z = 1 criticality, in g? = 2 spatial 
dimensions, it must have a linear-T dependence with a co-efficient that is a universal amplitude 
^4 X/ divided by the square of a non-universal velocity. 

er 

We will set c = 1 in the following. In general, we can think of the universal amplitude as a 
universal function A Xf = A Xf (Nf, N b ) for arbitrary Nf, N b . In the following we will calculate 
this functional dependence at large Nf, N b but for arbitrary values of the ratio N b /Nf. 

A Xf at N = oo: At iV = oo, we have a free theory of Dirac fermions and \f can be computed 
by simply evaluating the diagram in Fig.[8£a). 

2; 



Xj = -(Tr(n 2 )T^ [ ^Tr(-H + -k) 



We will evaluate the momentum integration first, using dimensional regularization, and then the 
frequency summation, using zeta function regularization. It is easy to obtain the final answer by 
evaluating the frequency summation before the momentum integration, in which case no cutoff or 
regularization is needed. However, on including gauge fluctuations, it appears more convenient to 
perform the momentum integration first, this N = oo calculation then provides an opportunity to 
illustrate the method by which we regulate the frequency sums later on. Performing the momentum 
integration in d spatial dimensions, we obtain 



We evaluate the frequency summation using 



9( s ) = T Y,r L Ts= 1T X - S ir- S (1 - 2- s ) C{s) (C4) 

From this, note that 

T]Tl = 0(0) = 
TyinlwJ = - % =Tln2 (C5) 
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Using these results in Eq. (IC3I) in taking the limit d — > 2, we obtain 

Xf = (Tr(T a ) 2 ) (C6) 

7T 

^4 X/ a? next to leading order in 1/N: We now turn to a computation of the diagrams of 
Fig. [8tb,c), these are the only contributions to at next to leading order. As will be clear be- 
low a full evaluation of these diagrams is not possible analytically. However we will be able 
to achieve our main goal, which is to extract the numerical value of the 1/N corrections to the 
N = oo value of A Xf computed above. 

A direct evaluation of the two diagrams, Fig. [8tb,c) can be written in the form (including the 
N = oo result above), 



Xf 



(Tr(r a ) 2 ) ^(0,0) + ^^ J 



d 2 k 
4^ 



D 3 (k,e r 



D 4 (k,e r 



D^k, e n ) D 2 (k, e n ) + {e 2 Jk 2 )D 1 {k, e n 

/d 2 a 
^-§Tr [2(4 (q, cj„)G*(q + k, u n + e n ) + G% (q, u n )G% (q + k, u n + e n )] 

"in 



(q, LU n )aiGl(q + k,u n + € n )a ] (<% - 



k 2 



where the functions have been introduced (Section Til Bl) and computed (Appendix lAl) before. 
Computations of £) 3 4 are presented in Appendix [Ql We note here that the dependence of A Xf on 
Nb enters through the iV& dependence of D\ 2 . A full numerical evaluation of A Xf was carried out 
and was found to reproduce the results of Sec. |IV]up to three significant digits. 



APPENDIX D: EVALUATION OF D 3 AND D 4 

In this Appendix we summarize our strategy to evaluate the functions D 3 and D 4 defined in 
Eq. (IC7I) . The techniques are similar to those used to evaluate D\ 2 and we will simply quote the 
final answers, 
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D 3 (k,e n ) = NTj^l 0Tr [2 (-iu; n + «t • q) 3 (-i(u n + e n ) + a ■ (q + k)) 1 

+ (-iu n -\- o • q) 2 (-i(a> n + e n ) + <7 • (q + k))~ 2 ] 
= / ctej^-2(:r - - (x - l)(5xe n + e n )u 3 n 

-(x-1) (2x 2 k 2 - xk 2 - 2x 2 e 2 n + 7xe 2 ) u 2 n 

— (x — 1) {7k 2 e n x 3 — 3e 3 l x 2 — 8k 2 e n x 2 + 3e 3 n x + 3k 2 e n x) u n 

-(x-1) ((x - l)x 2 k 4 + 3x 3 e 2 n k 2 - 2x 2 e 2 n k 2 - x 2 e 4 n ) 

((l-x)uj 2 n + x(uj n + e n ) 2 + ^x(l-x)f V ' 

This last lengthy expression was generated in ETgX by Mathematica, and also translated to Fortran 
by Mathematica. The integral over x was evaluated numerically, and the frequency summation was 
performed using the following expansion at large uo n 

1 k 2 - 4 -k A + 7k 2 e 2 n - 2e 4 n 

Similarly for _D 4 , we have 

D,(k, e n ) = U 3 - ^ W W 0Tr [2 (-iu n + a ■ q)' 3 a, (-i(u n + e n ) + a • (q + k))" 1 <r, 

+ (-iw n + cr • q)~ 2 Oi (-i(u n + e n ) + a ■ (q + k))~ 2 Gj] 
I diV -(1 - 3x)(x - l)u^ — (x - 1) (-10e n a; 2 + 3e n x + e n ) cj 3 



NT 



2tt 



- (x - 1) (-6A; 2 :r 3 - 8e 2 :r 3 - k 2 x 2 - 7e 2 n x 2 + 8k 2 x + 8e 2 n x) oo 2 n 



OJr. 



-(x-1) (-6k 2 e n x 4 - 8e 3 n x 3 + k 2 e n x 3 + e 3 n x 2 + 4k 2 e n x 2 + 3e 3 n x + 3k 2 e n x) 
-(x-1) (x 2 (x 3 - Ax 2 + Q x -3)k A - 3x 4 e 2 n k 2 + 6x 3 e 2 n k 2 - 2x 2 e 2 n k 2 - 2x 3 e 4 n + x 2 e 4 n ) 

((l- x )ul + x(u n + e n ) 2 + k 2 x(l-x)y K ) 

The frequency summand has the following expansion at large uo n 

-2A: 2 + e 2 n 4A; 4 - 17k 2 e 2 n - ej 
2ui + 6^ 
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